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La presente invention concerne un procede et un dispositif de pilotage 
d'un ensemble de restitution d'un champ acoustique comportant une piuralite 
d'elements de restitution, a partir d'une piuralite de signaux acoustiques ou au- 
diophoniques associes chacun a une direction generate de restitution predeter- 
minee, definie par rapport a un point de I'espace donne. 

Un tel ensemble de signaux est couramment designe par I'expression 
« signal multicanal » et correspond a une piuralite de signaux, appeles canaux, 
transmis en parallele ou multiplexes les uns avec les autres, chacun destine a un 
element ou un groupe d'elements de restitution, dispose dans une direction gene- 
rale predefinie par rapport a un point donne. 

Par exemple, un systeme multicanal classique et connu sous le nom 
de « 5.1 ITU-R BF 775-1 » et comporte cinq canaux destines a des elements de 
restitutions places dans cinq directions generates predetermines par rapport a 
un centre d'ecoute, definies par les angles 0°, + 30°, -30°, +1 10° et -110°. 

Une telle disposition correspond done a la disposition d'un haut-parleur 
ou un groupe de haut-parleurs devant au centre, un de chaque c6te devant a 
droite et a gauche et un de chaque c6te derriere a droite et a gauche. 

Les signaux de pilotage etant associes chacun a une direction deter- 
mine, ('application de ces signaux a un ensemble de restitution dont les ele- 
ments ne repondent pas a la configuration spatiale predetermine^ entraTne des 
deformations importantes du champ acoustique restitue. 

II existe des systemes qui integrant des moyens de retard sur les ca- 
naux, afin de compenser au moins partiellement, I'eloignement des elements de 
restitution par rapport au centre d'ecoute. Ces systemes ne permettent cepen- 
dant pas de prendre en compte la disposition dans I'espace de I'ensembie de 
restitution. 

II apparaTt done qu'aucun procede ou systeme existant, ne permet une 
restitution de bonne qualite a partir d'un signal de type multicanal avec un en- 
semble de restitution de configuration spatiale quelconque. 

a presen t e ^n vent iQ^^^o^^^Q^t-€te^ e n d 6 _ 
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finissant un procede et un systeme de pilotage de I'ensembie de restitution dont 
-la-confiymaUon spatiHle~6st que l conque 



L'invention a pour objet un procede de pilotage d'un ensemble de res- 
titution d'un champ acoustique comportant une piuralite d'elements de restitution 




associes chacun a une direction generate de restitution predeterminee definie par 
rapport a un point donne, pour obtenir un champ acoustique restitue de caracte- 
ristiques specifiques sensiblement independantes des caracteristiques intrinse- 
ques de restitution dudit ensemble, caracterise en ce qu'il comporte : 
5 - une etape de determination de caracteristiques au moins spatiales 

dudit ensemble de restitution, permettant la determination de parametres repre- 
sentatifs pour au moins un element dudit ensemble de restitution de sa position 
dans les trois dimensions de Tespace par rapport audit point donne ; 

- une etape de determination de filtres d'adaptation & partir desdites 
10 caracteristiques au moins spatiales dudit ensemble de restitution et desdites di- 
rections generates de restitution predetermine associees a ladite pluralite de 
signaux d'entree deformations acoustiques ; 

- une etape de determination d'au moins un signal de pilotage desdits 
elements dudit ensemble de restitution par I'application desdits filtres 

15 d'adaptation a ladite pluralite de signaux d'entree d'informations acoustiques; et 

- une etape de delivrance dudit au moins un signal de pilotage en vue • 
d'une application auxdits elements de restitution, 

Suivant d'autres caracteristiques : - 

- ladite etape de determination de caracteristiques au moins spatiales;,? 
20 dudit ensemble de restitution, comporte une sous-etape de saisie permettant de 

determiner tout ou partie des caracteristiques dudit ensemble de restitution ; 

- ladite etape de determination de caracteristiques au moins spatiales 
dudit ensemble de restitution, comporte une 6tape de calibrage permettant de 
delivrer tout ou partie des caracteristiques dudit ensemble de restitution ; 

25 - ladite sous-etape de calibrage comporte pour au moins un des ele- 

ments de restitution : 

- une sous-etape d'emission d'un signal specifique vers ledit au 
moins un element dudit ensemble de restitution ; 

- une sous-etape d'acquisition de Tonde sonore emise en re- 
30 ponse par ledit au moins un element ; 

- une sous-etape de transformation desdits signaux acquis en un 
nombre fini desdits coefficients representatifs de Tonde sonore emise ; et 
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- une sous-etape de determination de parametres spatiaux et/ou 
acoustiques dudit element a partir desdits coefficients representatifs de I'onde 
sonore emise ; 

- ladite sous-etape de calibrage comporte en outre une sous-etape de 
determination de la position dans au moins Tune des trois dimensions de i'espace 
dudit au moins un element dudit ensemble de restitution ; 

- ladite etape de calibrage comporte une sous-etape de determination 
de la reponse en frequence dudit au moins un element dudit ensemble de restitu- 
tion ; 

- ladite etape de determination de filtres d'adaptation comprend : 

- une sous-etape de determination d'une matrice de decodage 
representative de filtres permettant la compensation des alterations de restitution 
dues aux caracteristiques spatiales dudit ensemble de restitution ; 

- une sous-etape de determination d'une matrice de rayonnement 
multicanal ideale representative des directions generates predeterminees asso- 
cies a chaque signal d'informations de la pluralite des signaux d'entree ; et 

- une sous-etape de determination d'une matrice representative 
desdits filtres d'adaptation a partir de ladite matrice de decodage et de ladite ma- 
trice de rayonnement multicanal ; 

- ladite etape de determination de filtres d'adaptation comporte une 
pluralite de sous-etapes de calcul permettant de deiivrer un ordre limite de preci- 
sion spatiale des filtres d'adaptation, une matrice correspondant a une fenetre 
spatiale representative de la repartition dans I'espace de la precision souhaitee 
lors de la reconstruction du champ acoustique et une matrice representative du 
rayonnement de I'ensemble de restitution, ladite sous-etape de calcul de la ma- 
trice de decodage etant realisee a partir des resultats de ces sous-etapes de cal- 
cul ; 

- les matrices de decodage, de rayonnement multicanal ideal et 
d'adaptation sont independantes de la frequence, I'etape de determination d'au 

^oir4s^n^gnal^e^otage^sd4ts^efflents-dwlit ono o mb le-de-fes 
I'appiication desdits filtres d'adaptation correspondant a de simples combinaisons 
Hheairessuivies de retarcT ~ 



- ladite etape de determination de caracteristiques dudit ensemble de 
restitution permet la determination de caracteristiques acoustiques dudit ensem- 




ble de restitution et ledit procede comporte une etape de determination de fiitres 
de compensation de ces caract6ristiques acoustiques, ladite etape de determina- 
tion d'au moins un signal de pilotage comprenant alors une sous-etape 
d'application desdits fiitres de compensation acoustique ; 
5 - ladite etape de determination de caracteristiques acoustiques est 

adaptee pour delivrer des parametres representatifs pour au moins un element 
de sa reponse en frequence ; 

- ladite etape de determination d'au moins un signal de pilotage com- 
porte une sous-etape d'ajustement de gain et duplication de retards afin 

10 d'aligner temporellement le front d'onde des elements de restitution en fonction 
de leur distance par rapport audit point donn6. 

Uinvention a egalement pour objet un programme d'ordinateur com- 
prenant des instructions de code de programme pour Texecution des etapes du 
procede lorsque ledit programme est execute par un ordinateur. 

15 .Uinvention a encore egalement pour objet un support amovible du , 

type comprenant au moins un processeur de traitement et un element de me Tl , 
moire non volatile, caracterise en ce que ladite memoire comprend un pro-^ 
gramme comprenant des instructions de code pour Texecution des etapes du«. 
procede, lorsque ledit processeur execute ledit programme. 

20 Uinvention a encore pour objet un dispositif de pilotage d'un ensemble 

de restitution d'un champ acoustique comportant une pluralite d'elements de res- 
titution, comportant des moyens d'entree d'une pluralite de signaux d'entree 
d'informations acoustiques associes chacun a une direction g^nerale de restitu- 
tion predeterminee definie par rapport a un point donne, caracteris6 en ce qu'il 

25 comporte en outre : 

- des moyens de determination de caracteristiques au moins spatiales 
dudit ensemble de restitution, permettant la determination de parametres repre- 
sentatifs pour au moins un element dudit ensemble de restitution de sa position 
dans les trois dimensions de Tespace par rapport audit point donne ; 

30 - des moyens de determination de fiitres d'adaptation a partir desdites 

caracteristiques au moins spatiales dudit ensemble de restitution et des direc- 
tions generales de restitution pred6termin6e associees a ladite pluralite des si- 
gnaux d'entree d'informations acoustiques ; et 




- des moyens de determination d'au moins un signal de pilotage des- 
dits elements dudit ensemble de restitution par I'application desdits filtres 
d'adaptation a ladite pluralite de signaux d'entree d'informations acoustiques. 

Suivant d'autres caracteristiques de ce dispositif : 
5 - lesdits moyens de determination des caracteristiques au moins spa- 

tiales dudit ensemble de restitution comportent des moyens de saisie directe 
desdites caracteristiques ; 

- il est adapte pour §tre associe a des moyens de calibrage permettant 
la determination des caracteristiques au moins spatiales dudit ensemble de 

10 restitution ; 

- lesdits moyens de calibrage comprennent des moyens d'acquisition 
d'une onde sonore comportant quatre capteurs de pression disposes selon une 
forme generate de tetraedre ; 

- lesdits moyens de determination de caracteristiques sont adaptes 
1 5 pour la determination de caracteristiques acoustiques d'au moins un desdits ele- 
ments de restitution dudit ensemble de restitution, ledit dispositif comportant des 
moyens de determination de filtres de compensation acoustique a partir desdites 
caracteristiques acoustiques et lesdits moyens de determination d'au moins un 
signal de pilotage etant adaptes pour I'application desdits filtres de compensa- 

20 tion acoustique ; 

- lesdits moyens de determination des caracteristiques acoustiques 
sont adaptes pour la determination de la reponse en frequence desdits elements 
de I'ensemble de restitution. 

L'invention a egalement pour objet un appareil de traitement de don- 
25 nees audio et video comportant des moyens de determination d'une pluralite de 
signaux d'entree d'informations acoustiques associes chacun a une direction ge- 
nerate de restitution predetermine definie par un point donne, caracterise en ce 
qu'il comporte en outre un dispositif de pilotage d'un ensemble de restitution ; 

- lesdits moyens de determination d'une pluralite de signaux d'entree 
30 sont form es-dUme^ite^JeOecfare^^^ 

num§riques. 

CTrivention sera mieux^bmprise a la lectufe~de la description "aJFva"' 
suivre, donnee uniquement a titre d'exemple et faite en se referant aux dessins 
annexes, sur lesquels : 




- la Fig.1 est une representation d'un repere spherique ; 

- la Fig.2 est un schema cTun systeme de restitution selon ('invention ; 

- la Fig.3 est un organigramme du precede de I'invention ; 

- la Fig.4 est un schema de moyens de calibrage mis en oeuvre dans 
5 le procede de Tinvention ; 

- la Fig. 5 est un organigramme detaille de Tetape de calibrage ; 

- la Fig.6 est une representation simplifiee d'un capteur utilise pour la 
mise en oeuvre de Petape de calibrage ; 

- la Fig.7 est un organigramme detaille de I'etape de determination de 
10 filtres d'adaptation ; et 

- les Figs.8 et 9 sont des schemas de moyens de determination de sh 
gnaux de pilotage ; et 

- la Fig. 10 est un schema d'un mode de realisation d'un dispositif met- 
tant en oeuvre le procede de invention. 

15 Sur la figure 1, on a represents un repere spherique classique, de 

maniere a preciser le systeme de coordonnees auquel il est fait reference 
dans le texte. . F 

Ce repere est un repere orthonormal, d'origine O et comportant trpis 
axes {OX), {OY) et (OZ). .... 

20 Dans ce repere, une position notee x est decrite au moyen de ses 

coordonnees spheriques (r,0 9 $, oli r designe la distance par rapport a Torigine 0, 
6 I'orientation dans Je plan vertical et ^('orientation dans le plan horizontal. 

Dans un tel repere, un champ acoustique est connu si Ton d6finit en 
tout point S chaque instant t la pression acoustique not6e p(r,0 t $t), dont la trans- 
25 formee de Fourier temporelle est notee P(r, & t fif) oO/designe la frequence. 

L'invention est fondee sur ('utilisation d'une famille de fonctions spatio- 
temporelles permettant de decrire les caracteristiques de tout champ acoustique. 

Dans le mode de realisation decrit, ces fonctions sont des fonctions di- 
tes de Fourier-Bessel spheriques de premiere espece appelees par la suite fonc- 
30 tions de Fourier-Bessel. 

Dans une zone vide de sources sonores et vide d'obstacles, les fonc- 
tions de Fourier-Bessel sont solutions de ['equation des ondes et constituent une 
base qui engendre tous les champs acoustiques produits par des sources sono- 
res situees a I'exterieur de cette zone. 
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Tout champ acoustique tridimensionnel s'exprime done par une com- 
binaison lineaire des fonctions de Fourier-Bessel, selon I'expression de la trans- 
formee de Fourier-Bessel inverse qui s'exprime : 

5 Dar| s cette equation, les termes P it Jj) sont, par definition, les coeffi- 

cients de Fourier-Bessel du champ p(r,0,<f>,t), k=^-,c est la celerite du son dans 
I'air (340 ms' 1 ), ji(b-) est la fonction de Bessel spherique de premiere espece 
d'ordre / definie par j, (*)=^ J M2 (x) ou J v (x) est la fonction de Bessel de pre- 
miere espece d'ordre v, et y t m (0,fl est I'harmonique spherique reelle d'ordre / et 
1 0 de terme m, avec m allant de -/ a /, definie par : 



-4=~ ^''"'(cos 0) cos( m <p) pour m>0 

—^P t \zo%0) pour m=0 

■\I27T 

-4= P^'Ccos 0) sin( m tf>) pour m <0 



Dans cette equation, les P/ n (x) sont les fonctions de Legendre asso- 
ciees definies par : 



(I+m)l K X} dx"'* 1 W 



1 5 avec P,(x) les polynomes de Legendre, definis par : 

Les coefficients de Fourier-Bessel s'expriment aussi dans le domaine 
temporel par les coefficients Pl>m (t) correspondant a la transformee de Fourier 
temporelle inverse des coefficients Pijj). 
20 E n variante, le procede de I'invention utilise des bases de fonctions 
s'exprimant comme des combinaisons fineaires, eventuellement infinies, de fonc- 
^ions-d^Fotmer-Besseh — 



Sur la figure 2, on a represents schematiquement un systeme de resti- 



tution dans lequel on met^^uvrelen^ocede de I'invention. ~~~ 
25 c © systeme comporte un decodeur ou adaptateur 1 pilotant un en- 

semble de restitution 2 qui comporte une pluralite d'elements 3, a 3 Wl tels que 




des haut-parleurs, des enceintes acoustiques ou toute autre source sonore ou 
groupe de sources sonores, agences de maniere quelconque dans un lieu 
d'ecoute 4. On place arbitrairement, dans le lieu d'ecoute 4, Torigine O du repere 
que Ton appelle centre 5 de ['ensemble de restitution. 
5 L'ensemble des caracteristiques spatiales, acoustiques et electrody- 

namiques est considere comme les caracteristiques intrinseques de Tensemble 
de restitution 2. 

L'adaptateur 1 report en entree un signal SI de type multicanal compor- 
tant des informations acoustiques a restituer et un signal de definition SL compor- 
10 tant des informations representatives de caracteristiques au moins spatiales de 
Tensemble de restitution 2 et notamment permettant la determination de parame- 
tres representatifs pour au moins un element 3 n de Tensemble de restitution 2 de 
sa position dans les trois dimensions de Tespace par rapport au point donn6 5. 

A Tissue du traitement correspondant au procede de invention, 
1 5 l'adaptateur 1 emet & T attention de chacun des elements ou groupes d'elements 
3i a 3wde Tensemble de restitution 2, un signal sci a sc^de pilotage specifique. 

Sur la figure 3, on a represents schematiquement les etapes principa- 
les du procede selon Tinvention mis en ceuvre avec un systeme de restitution tel & 
que celui-decrit en reference a la figure 2. 
20 Ce procede comporte une etape 10 de determination de parametres 

de fonctionnement, adaptee pour permettre au moins la determination des carac- 
teristiques spatiales de ['ensemble de restitution 2. 

L'6tape 10 comporte une etape 20 de saisie des parametres et/ou une 
etape 30 de calibrage permettant de determiner et/ou mesurer des caracteristi- 
25 ques de Tensemble de restitution 2. 

Dans le mode de realisation d6crit, Tetape 10 comporte egalement une 
etape 40 de determination de parametres de description des directions generales 
predetermines associees aux differents canaux du signal d'entr£e multicanal SI. 

A Tissue de Tetape 10, des informations relatives au moins aux diffe- 
30 rentes directions generales predetermines associees a chacun des canaux 
d'entree ainsi qu'a la position dans les trois dimensions de Tespace de chacun 
des elements ou groupes d'elements 3 n de Tensemble de restitution 2, sont de- 
termines . 



• 




9 



10 



15 



20 



Ces informations sont utilises lors d'une etape 50 de determination 
des filtres d'adaptation permettant de prendre en compte les caracteristiques 
spatiales de I'ensemble de restitution 2 afin de definir des filtres d'adaptation du 
signal d'entree multicanal a la configuration spatiale specifique de I'ensemble de 
restitution 2. 

Avantageusement, I'etape 10 permet egalement de determiner des ca- 
racteristiques acoustiques pour tout ou partie des elements 3-, a 3 N de I'ensemble 
de restitution 2. 

Dans ce cas, le procede comporte une etape 60 de determination de 
filtres de compensation acoustique permettant de compenser I'influence des ca- 
racteristiques acoustiques specifiques des elements 3-i a 3 N . 

Les filtres definis lors des etapes 50 et avantageusement 60, peuvent 
ainsi §tre memorises, de sorte que les etapes 10, 50 et 60 ne doivent etre repe- 
tees qu'en cas de modification de la configuration spatiale de I'ensemble de resti- 
tution 2 et/ou de la nature du signal d'entree multicanal. 

Le procede comporte ensuite une etape 70 de determination des si- 
gnaux de pilotage sc-i a sc w destines aux elements de I'ensemble de restitution 2, 
comportant une sous-etape 80 d'application des filtres d'adaptation determines 
lors de I'etape 50 aux differents canaux c\(t) a c Q (t) formant le signal multicanal 
d'entree SI et avantageusement, une sous-etape 90 d'application des filtres de 
compensation acoustique determines lors de I'etape 60. 

Les signaux sc-i a sc N ainsi delivres, sont appliques aux elements 3-, a 
3 N de I'ensemble de restitution 2, afin de restituer le champ acoustique represen- 
ts par le signal multicanal d'entree SI avec une adaptation optimum aux caracte- 
ristiques spatiales et avantageusement acoustiques, de I'ensemble de restitution 
2. 

II apparaTt done que grace a la mise en ceuvre du procede de 
('invention, les caracteristiques du champ acoustique restitue sont sensiblement 
independantes des caracteristiques intrinseques de restitution de I'ensemble de 




2-^t-netamment^le-sa-configuration-spatiale: 

On va maintenant decrire plus en detail les etapes principales du 



procede de Pmvention. 



Lors de I'etape 20 de saisie de parametres un operateur ou un sys- 
teme de memoire adapte peut specifier tout ou partie des parametres de calcul et 
notamment : 

- des parametres x„ exprimes dans le repere spherique au moyen des 
5 coordonnees r n> 0* et fa , et representatifs de la position des elements 3„ par rap- 
port au centre d'ecoute 5 ; et/ou 

- des parametres H„(f), representatifs de la reponse en frequence des 
elements 3„ . 

Cette etape 20 est mise en ceuvre au moyen d'une interface de type 
10 classique telle qu'un micro-ordinateur ou tout autre moyen approprie. 

On va maintenant decrire plus en detail I'etape 30 de calibrage ainsi 
que des moyens de mise en ceuvre de cette etape. 

Sur la figure 4 on a represente le detail de moyens de calibrage. lis 
comportent un module 91 de decomposition, un module 92 de determination de 
15 reponse impulsionnelle et un module 93 de determination de parametres de cali7~ 
brage. ? 

Les moyens de calibrage sont adaptes pour etre connectes a un dis-;. 
positif d'acquisition sonore 100 tel qu'un microphone ou tout autre dispositif adap-i 
te, et pour etre connectes tour a tour a chaque element 3„ de I'ensemble de restiV? 
20 tution 2 afin de prelever des informations sur cet element. 

Sur la figure 5, on a represente le detail d'un mode de realisation de 
I'etape 30 de calibrage mise en oeuvre par les moyens de calibrage decrits pre- 
cedemment et permettant de mesurer des caracteristiques de I'ensemble de res- 
titution 2. 

25 Lors d'une sous-etape 32, les moyens de calibrage emettent un signal 

specifique u n (t) tel qu'une sequence pseudo-aleatoire MLS (Maximum Length 
Sequence) a I'attention d'un element 3„. Le dispositif d'acquisition 100 regoit, lors 
d'une sous-etape 34, 1'onde sonore emise par I'element 3„ en reponse a la recep- 
tion du signal u„(t) et transmet / signaux cp x (t) a cpfa) representatifs de I'onde re- 

30 gue au module de decomposition 91 . 

Lors d'une sous-etape 36, le module de decomposition 91 decompose 
les signaux captes par le dispositif d'acquisition 100 en un nombre fini de coeffi- 
cients de Fourier-Bessel qi, m (t). 
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Par exemple, le dispositif d'acquisition 100 est constitue de 4 capteurs 
de pression localises aux 4 sommets d'un tetraedre de rayon R ainsi que cela est 
represents en reference a la figure 6. Les signaux des 4 capteurs de pression 
sont alors notes cp,(f) a cp 4 (t). Les coefficients g 0 , 0 (t) a q hl (t) representatifs du 
5 champ acoustique capte se deduisent des signaux cp x (t) a cp 4 (i) selon les rela- 
tions suivantes : 

Q Qfi (f) = J_ cp M) + C W) + CPs(f) + CP 4 (f) 

3 

QlM = 8^? 2^jRf ( C *W + GP *W ~ CPs(f) - CP 4 (/)) 

Dans ces relations CPi(# a C7> 4 (/) sont les transformees de Fourier de 
10 c Pl (t) a cp 4 (t) et 2o,o(/) a £> 1(] tf) sont les transformees de Fourier de g Q>Q (t) a g u (t). 

Lorsque ces coefficients sont definis par le module 91, ils sont adres- 
ses au module 92 de determination de reponse. 

Lors d'une sous-etape 38, le module 92 de determination de reponse 
determine les reponses impulsionnelles h Plm (t) qui relient les coefficients de Fou- 
15 rier-Bessel g hm (t) et le signal emis u n (t). La methode de determination depend du 
signal specifique emis. Le mode de realisation decrit utilise une methode adaptee 
aux signaux de type MLS, comme par exemple la methode de correlation. 

La reponse impulsionnelle delivree par le module 92 de determination 
de reponse est adressee au module 93 de determination de parametres. 
20 Lors d 'une sous-etape 39, le module 93 deduit des informations sur 

des elements de I'ensemble de restitution. 

Dans le mode de realisation decrit, le module 93 de determination de 
parametres determine la distance /•„ entre I'element 3„ et le centre 5 a partir de sa 
reponse hp 0 ,o(t) et de la mesure du temps mis par le son pour se propager de 
25 I'element 3„ au dispositif d'acquisition 100, grace a des methodes d'estimation de 

retafd-su4Ma-repoRse-Ape^)(-4- ■ ■ 

La direction de I'element 3„ est deduite par calcul du maximum 
"o^Ta tran^i^e^eTbun^ 
hp\,\(t) prises a I'instant t oD hp 0 , 0 (t) presente un maximum. Avantageusement, les 
30 coordonnees 6„ et & sont estimees sur plusieurs instants, choisis de preference 




autour de I'instant ou &p 0 ,o(0 presente un maximum. La determination finale des 
coordonnees 0 n et est obtenue au moyen de techniques de moyennage entre 
ies differentes estimations. 

Ainsi, dans le mode de realisation decrit, le dispositif d'acquisition 100 
5 est en mesure d'encoder sans ambiguite I'orientation d'une source dans I'espace. 

En variante, les coordonnees 9 n et sont estimees a partir d'autres 
reponses parmi les hpiJf) dispbnibles ou sont estimees dans le domaine fre- 
quentiel a partir des reponses HPj tm (f), correspondent aux transformees de Fou- 
rier des reponses hpi tm (t). 
1 0 Ainsi Tetape 30 permet de determiner les parametres r„, B n et <f> iU 

Dans le mode de realisation decrit, le module 93 delivre 6galement la 
fonction de transfer! H n (f) de chaque el§ment 3„, a partir des reponses hpi, m (f) 
issues du module 92 de determination de reponse. 

Une premiere solution consiste & construire la reponse hp o,o(0 corres- 
15 pondant a la selection de la partie de la reponse hp 0t o(t) qui cpmporte un signal- 
non nul et denue des reflexions introduites par le lieu d'ecoute 4. La reponse en,'- 
frequence H n (j) est d6duite par transformee de Fourier de la reponse hp ' 0 ,o(0 pre-.* 
alablementfenetr^e. La fenetre peut §tre choisie parmi les fenetres classiques de % 
lissage, comme par exemple rectangulaire, Hamming, Hanning, et.Blackman. 
20 Une seconde solution, plus complexe consiste § appliquer un lissage 

sur le module et avantageusement sur la phase de la reponse en frequence 
HP 0 ,o(f) obtenue par transformee de Fourier de la reponse hp 0 ,o(t). Pour chaque 
frequence/, le lissage est obtenu par convolution de la reponse HP 0t o(f) par une 
fenetre centree sur / Cette convolution correspond a un moyennage de la re- 
25 ponse HP 0 ,o(f) autour de la frequence/. La fenetre peut etre choisie parmi ies fe- 
netres classiques, comme par exemple rectangulaires, triangles et Hamming. 
Avantageusement, la largeur de la fenetre varie avec la frequence. Par exemple, 
la largeur de la fendtre peut etre proportionnelle a la frequence / £ laquelle est 
appliquee le lissage. Par rapport a une fenetre fixe, une fenetre variable avec la 
30 frequence permet d'eiiminer au moins partiellement I'effet de salle dans Ies hau- 
tes frequences tout en evitant un effet de troncature de la reponse HP 0t o(f) dans 
les basses frequences. 

Les sous-6tapes 32 a 39 sont repetees pour tous les elements 3i a 3 W 
de Tensemble de restitution 2. 
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En variante, les moyens de calibrage comportent d'autres moyens 
d'acquisition d'informations relatives aux elements 3-, a Z N , tels que des moyens 
de mesure de position au laser, des moyens de traitement du signal mettant en 
ceuvre des techniques de formation de voies ou tout autre moyen approprie. 

Les moyens mettant en ceuvre l'etape 30 de calibrage sont constitues 
par exemple d'une carte electron ique ou d'un programme d'ordinateur ou de tout 
autre moyen approprie. 

L'etape 40 permet ainsi que cela a ete dit precedemment, de determi- 
ner des parametres decrivant le format du signal multicanal d'entree et notam- 
ment les directions predetermines generates associees a chaque canal. 

Cette etape 40 peut corresponds a une selection par un operateur 
d'un format parmi une liste de formats associes chacun a des parametres memo- 
rises, et peut egalement correspondre a une detection automatique de format 
effectuee sur le signal multicanal d'entree. Alternativement, le procede est adapte 
pour un unique format de signal multicanal donne. Dans encore un autre mode 
de realisation, l'etape 40 permet a un utilisateur de specifier son propre format en 
saisissant manuellement les parametres decrivant les directions associees a 
chaque canal. 

II apparatt que les etapes 20, 30 et 40 formant l'etape 10 de determi- 
nation de parametres, permettent au moins la determination de parametres de 
positionnement dans I'espace des elements 3„ de I'ensemble de restitution 2 et 
du format du signal multicanal SI. 

Sur la figure 7, on a represents un organigramme detaille de I'etape 50 
de determination des filtres d'adaptation. 

Cette etape comporte une pluralite de sous-etapes de calcul et de de- 
termination de matrices representatives des parametres determines prealable- 
ment. 

Ainsi, lors d'une sous-etape 51, un parametre L, appele ordre limite 
representatif de la precision spatiale souhaitee lors de l'etape 50 de determina- 
frdes4tttres^-adaptatiei^ 

- le plus petit angle a min forme par une paire d'elements de I'ensemble 
"deTestlution 2 est~calcuTS^5matiquement au moyen "d'une 7eTa¥oFtl^"onom§~ 
trique, telle que par exemple : 

a*\ »v = acosCsin^.siB^cosC^,,,-^) + cos^.cos^) 
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a min = min(a 7I / |7t 2) 

parmi I'ensemble des couples (nl, n2) tels que nl* n2\ et 
- ensuite, I'ordre maximum L est determine automatiquement comme 
etant !e plus grand entier respectant la relation suivant : 

L <n/ a min . 

L'etape 50 de determination de filtres d'adaptation comprend ensuite 
une sous-etape 52 de determination d'une matrice W^de ponderation du champ 
acoustique. Cette matrice correspond a une fendtre spatiale W{j-j) representa- 
tive de la repartition dans Tespace de la precision souhaitee lors de la reconstruc- 
tion du champ. Une telle fenetre permet de specifier la taille et la forme de la 
zone ou le champ doit etre correctement reconstruit. Par exemple, ii peut s'agir 
d'une boule centree sur le centre 5 de I'ensemble de restitution. Dans le mode de 
realisation d§crit, la fenetre spatiale et la matrice W sont independantes de la 
frequence. 

FT est une matrice diagonale de taille (ZH-1) 2 contenant des coefficients , 
de ponderation Wi et dans laquelle chaque coefficient W\ se trouve 2/+1 fois a la « 
suite sur la diagonale. La matrice Wa done la forme suivante : 

0 



w 0 0 



\ 0 

ow L 



Dans le mode de realisation decrit, les valeurs prises par les coeffi- 
cients Wi sont les valeurs d'une fonction telle qu'une fenetre de Hamming de taille 
de 2L+1 evailuee en /, de sorte que le param£tre W% est determine pour / allant de 
OaL 

L'etape 50 comporte ensuite une sous-etape 53 de determination 
d'une matrice M representative du rayonnement de I'ensemble de restitution no- 

tamment a partir des parametres de position x n . La matrice de rayonnement M 
permet de deduire des coefficients de Fourier-Bessel representant le champ 
acoustique qu'emet chaque element 3 n de I'ensemble de restitution en fonction 
du signal qu'il recoil 




M est une matrice de taille (L+l) 2 sur N, constitute d'elements Mj, m .„, 
les indices l,m designant la ligne l 2 +l+m et n designant la colonne n. La matrice M 
a done la forme suivante : 



0.0.1 
M 1,0,1 


M 0,0.2 

1 -A^ 1.-1,2*-- 
^1.0.2 

M ,,,,2 ••• 


•" -Wi.o.Ar 










M {..0.2'~ 








'■• M LtLtN 



Dans le mode de realisation decrit, les elements M,, m ,„ s'obtiennent a 
partir d'un modele de rayonnement en ondes planes, de sorte que : 

La matrice M ainsi definie est representative du rayonnement de 
I'ensemble de restitution. En particulier, M est representative de la configuration 
spatiale de I'ensemble de restitution. 

Les sous-etapes 51 a 53 peuvent etre executees sequentiellement ou 
simultanement. 

L'etape 50 de determination de filtres d'adaptation comprend ensuite 
une sous-etape 54 de prise en compte de I'ensemble des parametres du systeme 
de restitution 2 determines precedemment, afin de delivrer une matrice D de de- 
codage representative de filtres dits de reconstruction. 

En effet, les elements D nXm (f) de la matrice D correspondent a des fil- 
tres de reconstruction qui, appliques aux coefficients de Fourrier-Bessel P Km (f) 
d'un champ acoustique connu, permettent de determiner les signaux de pilotage 
d'un ensemble de restitution pour reproduire ce champ acoustique. 

La matrice de decodage D est done I'inverse de la matrice de rayon- 
nement M. 

La matrice D est obtenue a partir de la matrice Mau moyen de metho- 
de s d'inversion so us contrain t es faisant intervenir des parametres d'optimisation 
supplementaires. 

Dans le mode de realisation decrit, l'etape 50 est adaptee pour realiser 
une optimisation grace a la matrice de ponderation du champ acoustique FFqui 




permet notamment de reduire la distorsion spatiale dans le champ acoustique 
reproduit 

Cette matrice D est delivree notamment a partir de la matrice M, selon 
I'expression suivante : 
5 D = (M T (fM)" 1 mFw 

dans laquelle M T est la matrice transposee conjuguee de M. 

' Dans le mode de realisation dtcrit, les matrices Met Wsont indepen- 
dantes de la frequence, de sorte que la matrice D est egalement independante 
de la frequence. Elle est constitute d'elements notes D„j itn organises de la ma- 
10 niere suivante : 

A ao A t i,-i A, i t o Ajj * " ■ A,£ -l • ■ • A ,£.<>• " ■ A,x,l 
A,o,o A,i - 1 A,i ,o A, i,i • - * E>i,l-l * * • ■ • A,Z.,Z 

^AvAoA^lz-lAv^oA^U" * NyL—V 9 ^N y l t 0' ' % ^N t L,L^ 

Uetape 54 permet ainsi de delivrer la matrice Z> representative de fil- 
tres dits de reconstruction et permettant la reconstruction d'un champ acoustique 
a partir d'une configuration quelconque de Tensemble de restitution. Grdoe.d'* 

15 cette matrice, le procede de Tinvention permet de prendre en compte la configu- 
ration de Tensemble de restitution 2 et notamment de compenser les alterations ' 
du champ acoustique dues a sa configuration spatiale specifique. 

En variante, les parametres relatifs a ('ensemble de restitution 2 peu- 
vent dtre variables en fonction de la frequence. 

20 Par exempie, dans un tel mode de realisation, chaque element D n ^ m (f) 

de la matrice D peut etre determine en associant a chacun des iVsignaux de pilo- 
tage une fonction de directivite D n {9 9 <j>j) specifiant a chaque frequence / 
Tamplitude, et avantageusement la phase souhaitee sur le signal de pilotage sc n 
dans le cas d'une onde plane dans la direction 

25 On entend par fonction de directivite D n (0 9 &f) une fonction qui associe 

une valeur reelle ou complexe, eventuellement fonction de la frequence ou d'une 
plage de frequences, a chaque direction de Pespace. 

Dans le mode de realisation decrit, les fonctions de directivites sont in- 
dependantes de la frequence et notes D n (8,$. 

30 Ces fonctions de directivite D n (8,fy peuvent etre determinees en speci- 

fiant que certaines grandeurs physiques entre un champ ideal et le meme champ 
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reproduit par I'ensemble de restitution respectent des lois predeterminees. Par 
exemple, ces grandeurs peuvent etre la pression au centre et I'orientation du vec- 
teur Vitesse. Dans certains cas, on souhaite que seulement 3 signaux de pilotage 
soient actifs pour reproduire une onde plane. Les signaux de pilotage actifs, no- 
tes sc„i a sc n3 , sont ceux qui alimentent les elements de restitution dont les direc- 
tions sont les plus proches de la direction (0,$ de I'onde plane. Les elements de 
restitution actifs, notes 3„i a 3„ 3 , forment un triangle contenant la direction 
de I'onde plane. Dans ce cas, les valeurs des directivites D nJ (0,$ a D n3 (0,$ asso- 
cies aux 3 el6ments actifs 3 n1 a 3 n3 sont donnes par : 



10 q '~Ft=v 
avec 

(sin 6r> x cos tj> ni sin0„ 2 cos</>„ a sin& n3 cos0„A /sin0cos<^ 

sin0 ni sin0„, sin0„ 2 sm<j& n2 sin0 n3 sm<p„ 3 I r = I s in6>sin<j!> 

cosO ni cos0„ 2 cos 0 n3 j \ COS 9 

Dans cette relation, a correspond au vecteur contenant 
[D ni (e,®...D n3 (eM et les directions (4/, A/). (0^2) et (0„ 3 ,fa 3 ) correspondent 
1 5 respectivement aux directions des elements 3„i, 3„2 et 3„ 3 . 

On considere que les valeurs des directivites D n (0,® correspondant 
aux elements de restitution non actifs sont nulles. 

La relation precedente est repetee pour K directions d'ondes 
planes differentes. Ainsi, chacune des fonctions de directivite Z>„(0,$ est fournie 
20 sous la forme d'une liste de K echantilions. Chaque echantillon est fournit sous la 
forme d'un couple 011 (0 k ,$) est la direction de I'echantillon k 

et ou D n (0 k ,<fc) est la valeur de la fonction de directivite associee au signal de pilo- 
tage sc„ pour la direction 

Pour chaque frequence /, les coefficients D nJ Jf) de chaque fonction 
25 de directivite sont deduits des echantilions {((&,£), . Ces coefficients 
sont obtenus par inversion du processus d'echantillonnage angulaire qui permet 
de deduire les echantilions de la liste a partir d'une fonction 
de directivite fournie sous forme de coefficients d'harmoniques spheriq ues. Cette 
inversion peut prendre differentes formes afin de controler Interpolation entre ies 
-30- — eehantillons; - 

Dans d'autres modes de realisation, les fonctions de directivite sont di- 
rectement fournies sous la forme de coefficients D nJfm (f) de type Fourrier- Bessel. 
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Les coefficients D n j itn (f) ainsi determines sont utilises pour former la 

matrice 2>. 

L'etape 50 comporte ensuite une etape 55 de determination d'une ma- 
trice de rayonnement multicanale ideale S representative des directions genera- 
5 les predeterminees assoctees a chaque canal du signal multicanal d'entree SI. 

La matrice S est representative du rayonnement d'un ensemble de res- 
titution ideal, c'est a dire respectant parfaitement les directions generales prede- 
terminees du format multicanal. Chaque element Si tm . q (f) de la matrice S permet 
de deduire les coefficients de Fourier-Bessel Pi m (j) du champ acoustique ideale- 
1 0 ment restitue par chaque canal c q (t). 

La matrice S est d§terminee en associant a chaque canal d'entree c q (t) 
et avantageusement pour chaque frequence/, une figure de directivity represen- 
tative d'une distribution de sources supposees emettre le signal du canal c q (t). 

La distribution de sources est donnee sous la forme de coefficients 
15 d'harmoniques spheriques Si, m>q (j). Les coefficients Si^tf) sont ranges dans la 
matrice S de taille (£+l) 2 sur Q, oCi Q est le nombre de canaux. 

Dans le mode de realisation decrit, l'etape de mise en forme associe a. 
chaque canal c q (t) une source d'onde plane orientee dans la direction (6^).. 
correspondant a la direction (0j,$) associee au canal c q {t) dans le format multi- 
20 canal d'entree. Les coefficients Si initq (f) sont alors independants de la frequence, 
lis sont notes et s'obtiennent par la relation ; 

Dans d'autres modes de realisation, la matrice de rayonnement ideale 
S associe une distribution discrete de sources d'ondes planes a certains canaux 
25 pour simuler Peffet d'une ceinture de haut-parleurs. Dans ce cas, les coefficients 
Si, m , g sont obtenus par sommation des contributions de chacune des sources ele- 
mentaires. 

Dans encore d'autres modes de realisation, la matrice de rayonnement 
ideale S associe certains canaux c q (t) a Une distribution continue de sources 
30 d'ondes planes decrite par une fonction de directivite S q {9 9 <$. Dans ce cas, les 
coefficients Si >m>q de la matrice S sont obtenus directement par transformee de 
Fourier Spherique de la fonction de directivite S q (0,$. Dans ces modes de realisa- 
tion, la matrice S est independante de la frequence. 
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Dans d'autres modes de realisation plus complexes la matrice S asso- 
cie a certains canaux, une distribution de sources produisant un champ diffus. 
Dans ce cas, la matrice S varie avec la frequence. Ces modes de realisation sont 
adaptes aux formats multicanaux qui considered differemment les canaux avants 
et arrieres. Par exemple, dans les applications destinees a la restitution dans des 
salles cinema, les canaux arrieres sont souvent destines a recreer une ambiance 
diffuse. 

Dans d'autres modes de realisation, la matrice S associe a certains 
canaux des sources sonores dont la reponse n'est pas plate. Par exemple, dans 
le cas ou le format multicanal associe au canal c g (t) une source d'onde plane de 
reponse en frequence i/<% les S,,„, g (f) varient avec la frequence et s'obtiennent 
par la relation : 

Si le format multicanal associe a certains canaux une superposition 
des types de distributions de sources pre-citees, les coefficients S 1>m/1 (f) de la ma- 
trice de rayonnement s'obtiennent par sommation des coefficients associes a 
chaque type de distribution de source. 

Enfin, I'etape 50 comporte une sous-etape 56 de determination d'une 
matrice d'adaptation spatiale,* correspondant aux filtres d'adaptation a appiiquer 
au signal d'entree multicanal pour obtenir une restitution optimum prenant en 
compte la configuration spatiale de I'ensemble de restitution 2. 

La matrice d'adaptation spatiale A est obtenue a partir des matrices de 
mise en forme S et de decodage D au moyen de la relation : 

A=DS 

La matrice d'adaptation A permet de generer des signaux sa,(t) a sarffi 
adaptes a la configuration spatiale de I'ensemble de restitution a partir des ca- 
naux Cj (t) a c e (i). Chaque element A njq (j) est un filtre specifiant la contribution du 
canal c q {t) au signal adapte sa n (t). Grace a la matrice d'adaptation A, le precede 
de I'invention permet la restitution optimum du champ acoustique decrit par le 



signal multicanal' par un ensemble de restitution de configuration spatiale quel 
conque. 



Dans le mode de realisation decrit, les matrices D et S sont indepen- 
dantes de la frequence et la matrice A egalement. Dans ce cas, les elements de 
la matrice A sont des constantes notees A n ,„ et chacun des signaux adaptes sa,(t) 




a ttz*(0 s'obtient par de simples combinaisons lineaires des canaux d'entree c } (t) 
a c Q (t\ le cas echeant suivies de retard ainsi qu'il sera decrit ci-apres. 

Les filtres representes par la matrice A peuvent etre mis en oeuvre 
sous differentes formes de filtres et /ou methodes de filtrage. Dans le cas ou les 
5 filtres utilises sont parametres directement avec des reponses en frequence, les 
coefficients A n Jj) sont directement delivres par l'etape 50, Avantageusement, 
l'etape 50 de determination de filtres d'adaptation comprend une sous-etape 57 
de conversion afin de determiner les parametres des filtres pour d'autres metho- 
des de filtrage. 

10 Par exemple, les combinaisons de filtrage A n , q (f) sont converges en : 

- des reponses impulsionnelles finies a„, 7 (/) calcuiees par 
transformer de Fourier temporelle inverse de A ntg (f), chaque reponse inv 
pulsionnelle a ntg (t) est echantillonnee puis tronquee & une longueur propre 
a chaque reponse ; ou 
15 - des coefficients de filtres recursifs a reponses impulsiohnel- 

les infinies calcuiees a partir des A n>q (f) avec des methodes d'adaptation. rtl 
A Tissue de I'etape 50 les parametres des filtres d'adaptation A n>q (f) \, 
sont foumis. 

L'etape 60 permet ainsi que cela a ete dit precedemment, de determj-..:, 
20 ner les filtres de compensation des caracteristiques acoustiques des elements de 
I'ensemble de restitution 2 dans le cas ou des parametres relatifs a ces caracte- 
ristiques acoustiques tels que jes reponses en frequence H„(f), sont determines 
lors de l'etape 1 0 de determination des parametres. 

La determination de tels filtres, notes HtPtf) , a partir des reponses en 
25 frequence H„(f), peut §tre realisee de maniere classique en appliquant des me- 
thodes d'inversion de filtres, comme par exemple I'inversion directe, les metho- 
des de deconvolution, les methodes Wiener ou d'autres. 

En fonction des modes de realisation, la compensation porte unique- 
ment sur I'amplitude de la reponse ou encore sur I'amplitude et la phase. 
30 Cette etape 60 permet de determiner un filtre de compensation pour 

chaque element 3„ de I'ensemble de restitution 2 en fonction de ses caracteristi- 
ques acoustiques specifiques. 

De meme que precedemment, ces filtres peuvent §tre mis en ceuvre 
sous differentes formes de filtres et/ou methodes de filtrage. Dans le cas ou les 
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filtres utilises sont parametres directement avec des reponses en frequence, les 
reponses H$>{f) sont directement appliquees. Avantageusement, I'etape 60 de 
determination de filtres de compensation comprend une sous-etape de conver- 
sion afin de determiner les parametres des filtres pour d'autres methodes de fil- 
trage. 

Par exemple, les combinaisons de filtrage HfP(f) sont converties en : 

- des reponses impulsionnelles finies Af(0 calculees par 
transformee de Fourier temporelle inverse de HP(f) , chaque reponse im- 
pulsionnelle ftf>(t) est echantillonnee puis tronquee a une longueur propre 
a chaque reponse ; ou 

- des coefficients de filtres recursifs a reponses impulsionnel- 
les infinies calculees a partir des HP(f) avec des methodes d'adaptation. 

A Tissue de I'etape 60 les parametres des filtres de compensation 
HP(f) sontfournis. 

On va maintenant decrire plus en detail I'etape 70 de determination de 
signaux de pilotage. 

Cette etape 70 comporte une sous-etape 80 d'application des filtres 
d'adaptation represents par la matrices au signal multicanal d'entree SI corres- 
pondant au champ acoustique a restituer. Ainsi que cela a ete dit precedemment, 
les filtres d'adaptation A n ^(j) integrent les parametres caracteristiques de 
I'ensemble de restitution 2. 

Lore de la sous-etape 80, des signaux adaptes sa,(t) a sa^t) sont ob- 
tenus par Tapplication des filtres d'adaptation A„^(f) aux canaux cj(t) a c Q (t) du 
signal SI. 

Dans le mode de realisation decrit, la matrice d'adaptation A est inde- 
pendante de la frequence et les coefficients d'adaptation A„, g sont appliques de la 
facon suivante : 

Q 

%(4=5^(4^U - 

<7=3 

L'adaptation se pours uit par un ai usteme nt des gain^etj^pjlinaj^ 
de retards afin d'aligner temporellement les fronts d'onde des elements 3! a 3 W 
de I'ensemble de restitution 2 par rapport a ('element le plus eloigne. Les signaux 
adaptes sa t (i) a sa^t) se deduisent des signaux vj(*) a v^t) selon I'expression : 
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/*\ - , (* tnax( r n ) - n \ 
sa n {t) = n Vn\ t 



Dans d'autres modes de realisation, la matrice d'adaptation A varie 
avec la frequence et les filtres d'adaptation A„ jq (f) sont appliques de la fagon sui- 
vante : 

avec C q (f) la transform^ de Fourier temporelle du canal c q (t) et V n (f) 

defini par : 

y n (J) _ SAn (f) e -2*Jr*f/c 
Tn 

ou SAnif) est la transform^ de Fourier temporelle de sa n (i). 
10 Selon la forme des parametres des filtres d'adaptation A^qff), chaque 

filtrage des canaux c g (t) par les filtres d'adaptation A n>g (f) peut etre realise selon 
des methodes classiques de filtrage, telles que par exemple : 

- les parametres sont directement les reponses en frequence A nwq (J), et 
le. filtrage est effectue dans le domaine frequence!, par exemple, a I'aide des £ 

15 techniques usuelles de convolution par blocs ; 

- les parametres sont directement les reponses impulsionnelles finies,, 
<WX)i et le filtrage est effectue dans le domaine temporel par convolution ; ou 

- les parametres sont les coefficients de filtres recursifs a reponses 
impulsionnelles infinies, et le filtrage est effectue dans le domaine temporel au 

20 moyen des relations de recurrence. 

La sous etape 80 se termine par un ajustement des gains et 
Tapplication de retards afin d'aligner temporellement les fronts d'onde des ele- 
ments 3i a 3a/ de I'ensemble de restitution 2 par rapport a l'element le plus eloi- 
gne. Les signaux adaptes sa x (t) a sa^t) se deduisent des signaux v x (f) a v^{t) se- 

25 Ion ('expression : 

/A (. max( r n )-r n ^ 

SCn(t) = Vn V«l t V ^ I 

Sur la figure 8, on a represents la structure de filtrage correspondant a 
la sous etape 80 duplication des filtres d'adaptation spatiale telle que decrite 
precedemment 

30 Avantageusement, T6tape 70 comporte une sous-etape 90 de com- 

pensation des caracteristiques acoustiques de I'ensemble de restitution. Chaque 
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filtre de compensation HP(f) est applique au signal adapte sa„(f) correspondant 
afin d'obtenir le signal de pilotage sc n (t) de I'element 3„, selon la relation : 

SC n (f)=SA n (f) HPV) 

ou SC n (f) est la transformee de Fourier temporeile de sc„(t) et ou SA»(f) est la 
transformee de Fourier temporeile de sa n {t). 

L'application des filtres de compensation HP(f) des caracteristiques 
acoustiques est decrite en reference a la figure 9. 

Selon la forme des parametres de ces filtres, chaque filtrage des si- 
gnaux sa„(t) peut etre realise selon des methodes classiques de filtrage, telles 
que par exemple : 

- dans le cas ou les parametres de filtrage sont des reponses en fre- 
quence HP(f), le filtrage peut etre effectue au moyen de methodes de filtrage 

dans le domaine frequentiel, comme par exemple des techniques de convolution 
par bloc ; 

- dans le cas ou les parametres de filtrage sont des reponses impul- 
sionnelles h ( P(t), le filtrage peut etre effectue dans le domaine temporel par 
convolution temporeile ; 

- dans le cas 0C1 les parametres de filtrage sont des coefficients de re- 
lations de recurrence, le filtrage peut etre realise dans le domaine tempore! au 
moyen de filtres recursifs a reponse impulsionnelle infinie. 

Dans certains modes de realisation simplifies, le procede de I'invention 
ne compense pas les caracteristiques acoustiques specifiques des elements de 
Pensemble de restitution. Dans ce cas, I'etape 60 ainsi que la sous-etape 90 ne 
sont pas realisees et les signaux adaptes sai(t) a sa^t) correspondent directe- 
ment aux signaux de pilotage sci a sc N . 

Par l'application du procede de I'invention, chaque element 3-, a 3 N re- 
?oit done un signal de pilotage specifique sci a sc N et emet un champ acoustique 
qui contribue a la reconstruction optimale du champ acoustique a restituer. En 
i3ffetr4e-pilotage-sim 

construction optimale du champ acoustique correspondant au signal multicanal 
d'entree par re^sembleHjeliitiOT^ 

conque, soit encore ne correspond pas a une configuration fixe. 




Par ailleurs, d'autres modes de realisation du precede de ('invention 
peuvent etre envisages et notamment des modes de realisation inspires de tech- 
niques decrites dans la demande de brevet en France deposee le 28 fevrier 
2002, sous le n° 02 02 585. 
5 Notamment, Tetape 50 de determination des filtres d'adaptation spa- 

tiale peut prendre en compte de nombreux parametres d'optimisation tel que : 

- G n (f), representatif du gabarit de Telement 3„de Pensemble de restitu- 
tion specifiant la bande de frequence de fonctionnement de cet element ; 

- Ni.m.n(f)> representatif de la reponse spatio-temporelle de Telement 
10 3 n correspondant au champ acoustique produit dans le lieu d'ecoute 4 par 

Telement 3„, lorsque celui-ci regoit en entree un signal impulsionnel ; 

- W{r f j), decrivant pour chaque frequence / consideree une fenetre 
spatiale representative de la repartition dans Pespace de contraintes de recons- 
truction du champ acoustique, ces contraintes permettant de specifier la reparti- 

1 5 tion dans Tespace de Peffort de reconstruction du champ acoustique ; 

- Wi(j) t decrivant directement sous forme de ponderation des coeffi- t A 
cients de Fourier-Bessel et pour chaque fr6quence/consideree, une fenetre spa-^ 
tiale representative de la repartition dans I'espace de contraintes de reconstruc- 
tion du champ acoustique ; 

20 - R(f) t representatif, pour chaque frequence /consideree, du rayon de 

la fenetre spatiale lorsque celle-ci est une boule ; 

- /iff), representatif, pour chaque frequence /consideree, de la capaci- 
ty d'adaptation locale souhaitee a I'irregularite spatiale de la configuration de 
Pensemble de restitution ; 

25 - {(4, m k )}(f) 7 constituant pour chaque frequence/consideree, une liste 

de fonctions spatio-temporelles dont la reconstruction est imposee ; 

- L(f), imposant, pour chaque frequence / consideree, Pordre limite de 
determination de filtres ; 

- RM{f), definissant, pour chaque frequence / consideree, le modele de 
30 rayonnement des elements 3-j a Zn de Pensemble de restitution 2. 

Tout ou partie de ces parametres d'optimisation peuvent intervenir lors 
de la sous-etape 54 de determination de la matrice de decodage D. Ainsi, comme 
cela est decrit dans la demande de brevet en France deposee sous le numero 02 
02 585, les parametres Ni Mn (f) et RM(J) interviennent dans la sous-etape 53 de 
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determination de la matrice de rayonnement M, les parametres W(rJ), W x (f), R(f) 
interviennent dans la sous-etape 52 de determination de la matrice W, les para- 
metres {(h, m k )}(f) interviennent dans une sous-etape supplemental dans la 
determination d'une matrice F. La matrice de decodage D est alors determinee 
5 lors de la sous-etape 54, pour chaque frequence/, en fonction des matrices M, W 
et F et des parametres G„(f) et fx(J). 

Toujours, selon la demande de brevet 02 02 585, le calcul de la ma- 
trice D peut etre effectue frequence par frequence en considerant uniquement les 
elements actifs pour chaque frequence consideree. Cette methode de determina- 

10 tion de la matrice D fait intervenir le parametre G„(f) et permet d'exploiter au 
mieux un ensemble de restitution dont les elements ont des bandes de frequen- 
ces de fonctionnement differentes. 

II apparart que la mise en oeuvre du procede de I'invention decrit ici 
est plus efficace et done plus rapide que les procedes existants et notamment 

15 que le procede decrit dans la demande de brevet francais deposee sous le nu- 
mero 02 02 585. 

En effet, pour adapter un signal multicanal comportant Q canaux a un 
ensemble de restitution comportant N elements avec une precision spatiale 
d'ordre L, i\ apparatt que le procede de I'invention necessite QxN filtres 
20 d'adaptation au lieu des Q(L+l) 2 + (L+1) 2 N filtres necessaires a la mise en oeuvre 
du procede decrit dans la demande de brevet en France deposee sous le n° 02 
02 585. 

Par exemple, I'adaptation d'un signal « 5.1 ITU-R BF 775-1 » a un en- 
semble de restitution a 5 haut parleurs avec une precision d'ordre 5 necessite 25 
25 filtres au lieu de 360 filtres. 

Sur la figure 10, on a represente un schema d'un mode de realisation 
d'un appareil mettant en oeuvre le procede tel que decrit precedemment. 

Cet appareil comprend I'adaptateur 1 qui est forme d'une unite 110 de- 
livrant un signal multicanal telle qu'une unite de lecture de disques audio-video 
80— diteHecteur^VB-tt 2. Le s i gi l al ^tritirarrahdetiw 

elements de I'ensemble de restitution 2. Le format de ce signal SI est reconnu 
automatiquement paM^daplateLuM - ^ 

des parametres decrivant la direction generate predeterminee associee a chaque 
canal du signal SI. 




Selon I'invention, cet adaptateur 1 integre egalement une unite de cal- 
cul supplementaire 114 ainsi que des moyens de saisie d'informations 1 16. 

Par exemple, les moyens de saisie 116 sont formes d'une interface in- 
frarouge avec une telecommande ou encore avec un ordinateur et permettent a 
5 un utilisateur de determiner les parametres definissant les positions dans 
Pespace des elements de restitution 3i a 3^. 

Ces differents parametres sont utilises par le calculateur 1 14 pour de- 
terminer la matrice A definissant les filtres d'adaptation. 

Ulterieurement, le calculateur 114 applique ces filtres d'adaptation au 
10 signal multicanal SI afin de delivrer les signaux de pilotage sci a sc N a destination 
de I'ensemble de restitution 2. 

Bien entendu, le dispositif mettant en oeuvre I'invention peut prendre 
d'autres formes, telles qu'un logiciel mis en oeuvre sur un ordinateur ou encore 
un dispositif complet integrant des moyens de calibrage ainsi que des moyens de 
15 saisie et de determination des caracteristiques de I'ensemble de restitution plus 4 
complet. 

Ainsi, le procede peut egalement etre mis en oeuvre sous la forme % 
d'un dispositif dedie a ('optimisation de systemes de restitution multicanal, exte-.> 
rieur a un decodeur audio-video et associe a celui-ci. Dans ce cas, le dispositif 

20 est adapte pour recevoir en entree un signal multicanal et delivrer en sortie des 
signaux de pilotage d'elements d'un ensemble de restitution. 

Avantageusement, le dispositif est adapte pour etre connecte au dis- 
positif 100 d'acquisition necessaire a I'etape de calibrage et/ou est muni d'une 
interface permettant de saisir des parametres, notamment, la position des ele- 

25 ments de I'ensemble de restitution et eventuellement le format multicanal 
d'entn§e. 

Un tel dispositif d'acquisition 100 peut etre relie de maniere filaire ou 
non filaire (radio, infra-rouge) et peut etre integre a un accessoire, comme une 
telecommande, ou §tre independant. 
30 Le procede peut etre mis en oeuvre par un dispositif integre dans un 

element d'une chaTne audio-video charge du traitement de signaux multicanal, 
comme par exemple un processeur ou decodeur dit « surround », un amplifica- 
teur audio-video integrant des fonctions de decodage multicanal ou encore une 
chaTne audio-video compietement integree. 




Le precede de I'invention peut egalement etre mis en ceuvre dans une 
carte electronique ou dans une puce dediee. Avantageusement, il peut etre inte- 
gre sous la forme d'un programme dans un processeur de traitement de signal 
(DSP). 

5 Le procede peut prendre la forme d'un programme informatique desti- 

ne a etre execute par un ordinateur. Le programme regoit en entree un signal 
multicanai et delivre les signaux de pilotage d'un ensemble de restitution even- 
tuellement integre a cet ordinateur. 

Par ailleurs, les moyens de calibrage peuvent etre realises en mettant 
10 en ceuvre un procede different de celui decrit precedemment, tel que par exem- 
ple, un procede inspire de techniques decrites dans la demande de brevet en 
France deposee le 7 mai 2002 sous le numero 02 05 741 . 




REVEMD1CATIOMS 

1. Procede de pilotage d'un ensemble de restitution (2) d'un champ 
acoustique comportant une pluralite d'elements de restitution (3 n ) a partir d'une 
pluralite de signaux d'entree d'informations acoustiques (SI) associes chacun a 

5 une direction generate de restitution predeterminee d§finie par rapport a un point 
donne (5) de l'espace, pour obtenir un champ acoustique restitue de caracteristi- 
ques specifiques sensiblement independantes des caracteristiques intrinseques 
de restitution dudit ensemble (2), caracterise en ce qu'il comporte : 

- une etape (10) de determination de caracteristiques au moins spatia- 
10 les dudit ensemble de restitution (2), permettant la determination de parametres 

repr^sentatifs pour au moins un element (3 n ) dudit ensemble de restitution (2) de 
sa position dans les trois dimensions de l'espace par rapport audit point donne 
(5); 

- une etape (50) de determination de filtres d'adaptation (A) a partir 
15 desdites caracteristiques au moins spatiales dudit ensemble de restitution (2) et 

desdites directions generates de restitution predeterminee associees a ladite plu- 
ralite de signaux d'entr6e d'informations acoustiques (SI) ; 

- une etape (70) de determination d'au moins un signal de pilotage- 
desdits elements dudit ensemble de restitution par Implication desdits filtres 

20 d'adaptation a ladite pluralite de signaux d'entree d'informations acoustiques (SI); 
et 

- une etape de delivrance dudit au moins un signal de pilotage en vue 
d'une application auxdits Elements de restitution (3 n ). 

2. Procede selon la revendication 1 , caracterise en ce que ladite etape 
25 (10) de determination de caracteristiques au moins spatiales dudit ensemble de 

restitution (2), comporte une sous-etape (20) de saisie permettant de determiner 
tout ou partie des caracteristiques dudit ensemble de restitution (2). 

3. Procede selon Tune quelconque des revendications 1 et 2, caracte- 
rise en ce que ladite 6tape (10) de determination de caracteristiques au moins 

30 spatiales dudit ensemble de restitution (2), comporte une 6tape (30) de calibrage 
permettant de delivrer tout ou partie des caracteristiques dudit ensemble de 
restitution (2). 
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4. Precede selon la revendication 3, caracterise en ce que ladite sous- 
etape de calibrage (30) comporte pour au moins un des elements de restitution 

(3„): 

- une sous-etape d'emission (32) d'un signal specifique (u n (t)) vers le- 
dit au moins un element (3 n ) dudit ensemble de restitution (2) ; 

- une sous-etape d'acquisition (34) de I'onde sonore emise en reponse 
par ledit au moins un element (3 n ) ; 

- une sous-etape de transformation (36) desdits signaux acquis en un 
nombre fini desdits coefficients representatifs de I'onde sonore emise ; et 

- une sous-etape (39) de determination de parametres spatiaux et/ou 
acoustiques dudit element (3 n ) a partir desdits coefficients representatifs de 
I'onde sonore emise. 

5. Procede selon Tune quelconque des revendications 3 et 4, caracte- 
rise en ce que ladite sous-etape de calibrage (30) comporte en outre une sous- 
etape de determination de la position dans au moins Tune des trois dimensions 
de I'espace dudit au moins un element (3 n ) dudit ensemble de restitution (2). 

6. Procede selon I'une quelconque des revendications 3 a 5, caracteri- 
se en ce que ladite etape de calibrage (30) comporte une sous-etape de determi- 
nation de la reponse en frequence (H n (f)) dudit au moins un element (3 n ) dudit 
ensemble de restitution (2). 

7. Procede selon I'une quelconque des revendications 1 a 6, caracteri- 
se en ce que ladite etape (50) de determination de filtres d'adaptation comprend : 

- une sous-etape (54) de determination d'une matrice de decodage (D) 
representative de filtres permettant la compensation des alterations de restitution 
dues aux caracteristiques spatiales dudit ensemble de restitution (2) ; 

- une sous-etape (55) de determination d'une matrice de rayonnement 
multicanal ideale (S) representative des directions generates predeterminees as- 
socies a chaque signal d'informations de la pluralite des signaux d'entree (SI) ; et 

- une sous-etape (56) de determination d'une matrice (A) representa- 
twe-desdtts-filtres^adap^ 



ladite matrice de rayonnement multicanal (S). 

8. Procede seloFliTrevendTcatio^ 
(50) de determination de filtres d'adaptation comporte une pluralite de sous- 
etapes de calcul (51, 52, 53) permettant de delivrer un ordre (L) limite de la preci- 
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sion spatiale des filtres d'adaptation, une matrice (W) correspondant a une fene- 
tre spatiale representative de la repartition dans I'espace de la precision souhai- 
tee lors de la reconstruction du champ acoustique et une matrice (M) representa- 
tive du rayonnement de I'ensemble de restitution (2), ladite sous-etape (54) de 
5 calcu! de ia matrice (2) de ddcodage (D) etant realis6e a partir des resultats de 
ces sous-etapes de calcul. 

9. Procede selon Tune des revendications 7 ou 8, caracterise en ce 
que les matrices de decodage (D), de rayonnement multicanal ideal S et d'adap- 
tation A sont independantes de la frequence, I'dtape (70) de determination d'au 

10 moins un signal de pilotage desdits elements dudit ensemble de restitution par 
{'application desdits filtres d'adaptation correspondant a de simples combinaisons 
lineaires suivies de retard. 

10. Proced6 selon Tune quelconque des revendications 1 a 9, caracte- 
rise en ce que ladite etape (10) de determination de caracteristiques dudit en- 

1 5 semble de restitution (2) permet la determination de caracteristiques acoustiques 
dudit ensemble de restitution (2) et en ce que ledit procede comporte une etape 
(60) de determination de filtres de compensation de ces caracteristiques acousth 

ques, ladite etape (70) de determination d'au moins un signal de pilotage com- 

'*, 

prenant alors une sous-etape (90) d'appiication desdits filtres de compensation 
20 acoustique. 

11. Procede selon la revendication 10, caracterise en ce que ladite 
etape (10) de determination de caracteristiques acoustiques est adaptee pour 
delivrer des parametres representatifs pour au moins un element (3 n ) de sa re- 
ponse en frequence (H n (f)). 

25 12. Procede selon Tune quelconque des revendications 1 a 1 1 , carac- 

terise en ce que ladite etape (70) de determination d'au moins un signal de pilo- 
tage comporte une sous-etape d'ajustement de gain et d'appiication de retards 
afin d'aligner temporellement le front d'onde des elements de restitution (3 n ) en 
fonction de leur distance par rapport audit point donne (5). 

30 13. Programme d'ordinateur comprenant des instructions de code de 

programme pour I'execution des etapes du procede selon Tune quelconque des 
revendications 1 a 12 lorsque ledit programme est execute par un ordinateur. 

14. Support amovible du type comprenant au moins un processeur de 
traitement et un element de memoire non volatile, caracterise en ce que ladite 




memoire comprend un programme comprenant des instructions de code pour 
I'execution des etapes du procede selon Tune quelconque des revendications 1 a 
12, lorsque ledit processeur execute ledit programme. 

15. Dispositif de pilotage d'un ensemble de restitution (2) d'un champ 
5 acoustique comportant une plurality d'§lements de restitution (3 n ), comportant 

des moyens (112) d'entree d'une pluralite de signaux d'entree deformations 
acoustiques (SI) associes chacun & une direction generate de restitution prede- 
terminee definie par rapport a un point donne (5), caracterise en ce qu'il com- 
porte en outre : 

10 - des moyens (116) de determination de caracteristiques au moins 

spatiales dudit ensemble de restitution (2), permettant ia determination de para- 
metres representatifs pour au moins un element (3 n ) dudit ensemble de restitu- 
tion (2) de sa position dans les trois dimensions de I'espace par rapport audit 
point donne (5) ; 

15 - des moyens (114) de determination de filtres d'adaptation (A) a partir 

desdites caracteristiques au moins spatiales dudit ensemble de restitution (2) et 
des directions generates de restitution predeterminee associees a ladite pluralite 
des signaux d'entree d'informations acoustiques (SI) ; et 

- des moyens (1 14) de determination d'au moins un signal de pilotage 

20 (sc n ) desdits elements (3 n ) dudit ensemble de restitution (2) par I'application des- 
dits filtres d'adaptation (A) a ladite pluralite de signaux d'entree d'informations 
acoustiques (SI). 

16. Dispositif selon la revendication 15, caracterise en ce que lesdits 
moyens de determination des caracteristiques au moins spatiales dudit ensemble 

25 de restitution (2) comportent des moyens (116) de saisie directe desdites carac- 
teristiques. 

17. Dispositif selon Tune quelconque des revendications 15 et 16, ca- 
racterise en ce qu'il est adapte pour etre associe a des moyens de calibrage (91, 
92, 93, 100) permettant la determination des caracteristiques au moins spatiales 

-30 dudit-ensem bl e-de^esiitutian_£2) - — 

18. Dispositif selon la revendication 17, caracterise en ce que lesdits 
moyensHfle calibrage compi^Trerrird es ir iuy eris^'-acquisition-d^na-onde-sonere 
(100) comportant quatre capteurs de pression disposes seion une forme generate 
de tetraedre 



^ 

19. Dispositif selon Tune quelconque des revendications 15 a 18, ca- 
racterise en ce que lesdits moyens de determination de caracteristiques sont 
adaptes pour la determination de caracteristiques acoustiques d'au moins un 
desdits elements (3 n ) dudit ensemble de restitution (2), ledit dispositif comportant 

5 des moyens de determination de filtres de compensation acoustique a partir des- 
dites caracteristiques acoustiques et lesdits moyens de determination d'au moins 
un signal de pilotage etant adaptes pour ('application desdits filtres de compensa- 
tion acoustique. 

20. Dispositif selon la revendication 19, caracterise en ce que lesdits 
10 moyens de determination des caracteristiques acoustiques sont adaptes pour la 

determination de la reponse en frequence (H n (f)) desdits elements (3 n ) de 
I'ensemble de restitution (2). 

21. Appareil de traitement de donnees audio et video comportant des 
moyens (112) de determination d'une pluralite de signaux d'entree d 'informations 

15 acoustiques (SI) associes chacun a une direction g6nerale de restitution prede- 
terminee definie par un point donne (5), caracterise en ce qu'il comporte en outre 
un dispositif de pilotage d'un ensemble de restitution (2) selon Tune quelconque 
des revendications 1 a 19. 

22. Appareil selon la revendication 21, caracterise en ce que lesdits 
20 moyens de determination d'une pluralite de signaux d'entree sont formes d'une 

unite (112) de lecture et de decodage des disques audio et/ou video numeriques. 
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